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Abstract 

We propose a new test for the hypothesis that a bivariate copula is an Archimedean 
copula. The test statistic is based on a combination of two measures resulting from the 
characterization of Archimedean copulas by the property of associativity and by a strict 
upper bound on the diagonal by the Frechet-upper bound. We prove weak convergence of 
this statistic and show that the critical values of the corresponding test can be determined by 
the multiplier bootstrap method. The test is shown to be consistent against all departures 
from Archimedeanity if the copula satisfies weak smoothness assumptions. A simulation study 
is presented which illustrates the finite sample properties of the new test. 
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•Eh 1 Introduction 

X 

Let F be a bivariate continuous distribution function with marginal distribution functions F\ and 



F 2 . By Sklar's Theroem [see Sklarj (1959)] we can decompose F as follows 



F(x) =C(F 1 (x 1 ),F 2 {x 2 )), x=(x b x 2 )GM 2 , (1.1) 

where C is the unique copula associated to F. By definition, C is a bivariate distribution function 
on the unit square [0, l] 2 whose univariate marginals are standard uniform distributions on the 



interval [0,1]. Equation (1.1) is usually interpreted in the way that the copula C completely 
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characterizes the information about the stochastic dependence contained in F. For an extensive 



exposition on the theory of copulas we refer the reader to the monograph Nelsen (2006). 



In the last decades, various parametric models for copulas have been developed, among which 



the class of Archimedean copulas forms one the most famous and largest class, see Genest and 



MacKay 


(1986 


); 


Nelsen 


(2006 


); 



used copulas, such as Clayton-, Gumbel- and Frank-copulas are in fact Archimedean copulas. The 
elements of this class may be characterized by a continuous, strictly decreasing and convex function 
$ : [0, 1] -»■ [0, oo] satisfying $(1) = such that 

C(u) = $ M1 [$(«i) + $(u 2 )] for all u = (u u u 2 ) G [0, l] 2 . 

The function $ is called the generator of C and its pseudo-ineverse is defined as the 

usual inverse for t G [0, $(0)] and is set to for t > $(0). The prominence of the class of 

Archimedean copulas basically stems from the fact that they are easy to handle and to simulate, 
see Genest et al. (2011). While the estimation of Archimedean copulas has been investigated in 



Genest and Rivest (1993) and recently more thoroughly in Genest et al. (2011), the issue of testing 
for the hypothesis that the copula is an Archimedean one has found much less interest in the 
literature. The present paper fills this gap by developing a consistent test for this hypothesis. 



Our interest in this problem stems from recent work of Genest and Rivest (1993), Wang and Wells 



(2000) and Naifar ( 2011[ ) who proposed Archimedean copulas for modeling dependencies between 



bivariate observations (among many others). We also refer to the work of Rivest and Wells (2001) 
who used Archimedean copulas for modeling the dependence in the context of censored data. 



To the best of our knowledge, the only available test hitherto has been discussed in Jaworski (2010). 
This author proposed a procedure which is based on a characterization of Archimedean copulas 
similar to the one stated in Theorem 4.1.6 in Nelsen (2006) [which dates back to Ling (1965)]. To 



be precise recall that a bivariate copula C is called associative if and only if the identity 

C(x,C(y,z))=C(C(x,y),z) 



;i-2) 



holds for all (x,y,z) G [0, l] 3 . Theorem 4.1.6 in Nelsen (2006) shows that a bivariate copula C is 
an Archimedean copula if and only if C is associative and the inequality C(u,u) < u holds for 
all u G (0, 1), i.e. on the diagonal C is strictly dominated by the Frechet-upper bound M(u) = 
min(ui, U2). The procedure suggested in Jaworski (2010) is in fact to test for associativity in order 



to check the validity of an Archimedean copula model. The corresponding test statistic is defined 

as 



T n (x, y) = Vn(C n (x, C n (y, y)) - C n (C n (x, y),y), 

where (x,y) is some fixed point in the open cube (0, l) 2 and C n denotes the empirical copula, see 
Section [2] for details. The main advantage of this approach is its simplicity, in particular the simple 
limit distribution of the resulting test statistic, which is in fact normal. On the other hand this 



simplicity has its price in terms of consistency. In our opinion, the method proposed by Jaworski 



(2010) has at least three mayor drawbacks. First of all, it is clearly not consistent against a large 
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class of alternatives since it only tests for equation (1.2) with y = z. Second, Jaworski (2010) 
uses a pointwise approach in order to test for a global hypothesis as in (1.2). This means that the 
test may not reject the hypothesis because (1.2) is satisfied at the particular point (x,y,y) under 
investigation, although there may exist many other points where (1.2) is not satisfied. Third, there 
exist copulas which are in fact associative but not Archimedean. These problems also have strong 
implications for the practical applicability of the test as demonstrated by results in a simulation 
study in Jaworski (2010), where the sample size has to be chosen extremely large in order to get 
reasonable rejection probabilities. 

To the best of our knowledge there exists no test for an Archimedean copula, which is consistent 
against general alternatives and it is the primary purpose of this paper to develop such a proce- 
dure and to investigate its statistical properties. We propose a test statistic which is based on a 
combination of two measures resulting from the characterization of Archimedean copulas, namely 
the property of associativity as described in (1.2) and the strict upper bound on the diagonal 
C(u,u) < u for all u G (0, 1). In Section [2] we define a new process which is based on an estimate 



of the difference of the left and right hand side of the defining equation (1.2) for associativity. We 
prove weak convergence of this process in the space of all uniformly bounded functions on the cube 
[0, l] 3 . As a consequence, we also obtain weak convergence of a corresponding Cramer-von-Mises 
and a Kolmogorov-Smirnov type statistic. Because the asymptotic distribution depends in a com- 
plicated manner on the underlying copula we propose a multiplier bootstrap procedure to obtain 
the critical values and show its validity. As a first main result we obtain a test for associativity, 
which is consistent against all alternatives satisfying weak smoothness assumptions on C. In Sec- 
tion [3] we utilize these findings to develop an asymptotic test for the hypothesis of Archimedeanity. 
Finally in Section [4] we investigate the finite sample performance of the new test by means of a 
simulation study. 



2 Testing Associativity 

2.1 The test statistic and its asymptotic behavior 

In the following let Xi, . . . ,X n , Xj = (Xn, X i2 ) denote independent identically distributed bivari- 
ate random vectors with continuous distribution function F, marginal distribution functions F% 
and F 2 and copula C = F(F± , F 2 ~). In this paragraph we will introduce a test statistic for the 



null hypothesis that the underlying copula is associative, i.e. C satisfies condition (1.2) for all 
(x,y,z) G [0,1] 3 . 

For this purpose we briefly summarize relevant notations and results on the empirical copula, which 
is the simplest and most popular nonparametric estimator of the copula. In particular we define 
the empirical copula by 

C n (u) = F n (F- 1 (u 1 ),F- 2 (u 2 )), 

where F n (x) = n~ x XT=i < x} and F np (x p ) = n' 1 Y^=i ^-{X P < %p},P = 1,2 are the joint 
and marginal empirical distribution functions of the sample Xi, . . . , X n , respectively. It is a well 
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known result that under the assumptions of continuous partial derivatives of C the corresponding 
empirical copula process 



C, 



n(C n - C) 



(2.1) 



nian et al. 


(2004 


); 


Tsukahara 


(2005 



converges weakly towards a Gaussian limit field Gc in /°°([0, l] 2 ), see Riischendorf (1976); Ferma 

(2005) among others. Defining C p as the p-th partial derivative of 



(2.2) 



C (p — 1, 2) the process Gc can be expressed as 

G c (x) = B c (x) - C r i(x)l (ari, 1) - C' 2 (x)B c (l, x 2 ) 



with the copula-brownian bridge Be, i.e. Be is a centered Gaussian field with Cov(Bc(x),Bc(y)) = 
C(xAy) — C(x)C(y), where the minimum of two vectors is defined component-wise. As explained 
in Segers (2011) the assumption of continuity of the partial derivatives of C on the whole unit 
square does not hold for many (even most) commonly used copula models and as a consequence 
Segers provides the result that the following nonrestrictive smoothness condition is sufficient in 



order to obtain weak convergence of the empirical copula process defined in (2.1). 



Condition 2.1. For p = 1, 2 the first order partial derivative C p of the copula C with respect to 
x p exists and is continuous on the set V p = {u 6 [0, l] 2 : < u p < 1}. 

Now, in order to test for associativity we consider the process 

B. n (x,y, z) = y/n{C n (x,C n (y,z)) - C n (C n (x,y), z)} , 



where (x,y,z) G [0, l] 3 . The asymptotic properties of the process {H n (x, y, z)}( !B) j /j «)g[o,i]3 are sum- 
marized in the following Theorem. Throughout this paper l°°(T) denotes the set of all uniformly 
bounded functions on T, and the symbol denotes uniform convergence in a metric space (which 
will be specified in the corresponding statements). 



Theorem 2.2. If the copula C is associative and satisfies Condition 2A_, then it holds 

U n U c m r°([o,i] 3 ), 

where the limit field He can be expressed as 
Mc{x, y, z) = G c (x, C(y, z)) - G c (C(x, y),z) + C 2 (x, C(y, z))G c (y, z) - C 1 {C{x ) y), z)G(x, y). 

Proof. If the copula C is associative we can write the process M. n as 

U n = v^{<£(C n )-$(C)}, 
where the functional $ : D$ — > /°°([0, l] 3 ) is defined for 

a E D$ = {F : F cdf on [0, l] 2 } 
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by 

$(a)(x, y, z) = a(x, a(y, z)) - a(a(x, y),z). 



We will show later that under Condition 2.1 the mapping $ is Hadamard-differentiable at C 
tangentially to the space 

O = {7 G C[0, l] 2 I 7(u) = for all u G [0, l] 2 s.t. C(u) G {0, 1}} , 

with derivative given by 

& c (a)(x,y, z) = a(x,C(y, z)) - a(C(x,y),z) + C 2 (x, C(y, z))a(y, z) - C 1 (C(x,y), z)a(x,y). 



Observing that Be G D a.s., the functional delta method, see Theorem 3.9.4 in Van der Vaart 



and Wellner (1996), yields the assertion. 

We now briefly sketch how to see the Hadamard-differentiability of the mapping $: let t n — > and 
a n G Z°°([0, l] 2 ) with a n a G D such that C + t n a n G D$. Then 



^'{^(C + ^«n) - $(C)} = L nl + L n2 - L 
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where 

L n i(x,?/,z) = a n (x, (C + t„a n )(y,2)) - a n ((C + t n a n )(x, y), z) 
L n2 (x, y, z) = t~ l {C(x, (C + t n a n )(y, z)) - C(x, C(y, z))} 
L n3 (x, y, z) = t- l {C((C + t n a n )(x, y),z)- C(C(x, y), z)}. 

Exploiting the fact that a n converges uniformly to a bounded function and that a is uniformly 
continuous one can conclude that L n i(x,y,z) = a(x,C(y, z)) — a(C(x,y), z) + o(l) uniformly in 
(x,y,z) G [0, l] 3 . Regarding the summand L n2 we have to split the investigation in two cases. 
First, we consider all those (x,y,z) G [0, l] 3 for which C(y,z) G (0,1). A Taylor expansion of 
C(x, •) at C(y,z) yields 

L n2 (x, y, z) = C 2 (x, C(y, z))a n (y, z) + r n (x, y, z), 

where the error term can be written as 

r n (x,y,z) = (C 2 (x,u n ) - C 2 {x, C(y, z)))a n (y, z) 

with some intermediate point u n between C(y, z) and (C + t n a n )(y, z). The main term uniformly 
converges to C 2 (x,C(y, z))a(y, z) [note that partial derivatives of copulas are uniformly bounded 
by 1] and it remains to show that r n (x, y, z) = o(l) uniformly in (x, y, z) with C(y, z) G (0, 1). 
To see this, we will show at the end of this proof that for any e > there exists a 5 > 0, such that 

limsup sup |a n (v)| < e. (2.3) 
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where v = (y,z), As = {v G [0, 1] 2 1 C(v) G [0,5) U (1 — 8, 1]}. Then, since partial derivatives of 
copulas are bounded by 1, we can conclude that 

limsup sup \r n (x,y, z)\ < e. 

n->oo xe[0,l],(y,z)eA s 



Due to Condition 



2.1 



the partial derivative C% is uniformly continuous on the quadrangle [0, 1] x 
[5, 1 — 5}. Thus, since a is uniformly bounded and since u n — > C(y, z), we obtain uniform conver- 
gence of r n (x, y, z) to for all (y, z) s.t. C(y, z) G [5, 1 — 5], i.e. for (y, z) G [0, l] 2 \ Ag. Combining 
the two facts derived above, it follows that 

limsup sup \r n (x,y, z)\ < e. 

n-s>oo xe[0,l},C{y,z)e(0,l) 

Since e > was arbitrary, this limsup must be zero. Summarizing, the case (x, y, z) G [0, l] 3 such 
that C(y, z) G (0, 1) is finished. 

In the remaining case C(y, z) G {0, 1}, i.e. (y, z) G A , Lipschitz-continuity of C entails that 

\L n2 (x,y,z)\ =t; i 1 \C(x,C(y,z) +t n a n (y,z)) - C(x,C(y, z))\ < a n {y,z) = a(y,z) + o(l) = o(l) 



uniformly in (x,y,z) since in this case a(y,z) 
analogously. 



0. Finally, the summand L n3 may be treated 



To complete the proof it remains to show (2.3). Exploiting uniform convergence of a n , uniform 
continuity of a and the fact that a(v) = for all v G Aq = {v | C(v) G {0, 1}}, we can conclude 
that there exists a k > such that |a n (v)| < e for all v G Aq = {v | 3 u G Aq s.t. ||u — v|| < n} 
and sufficiently large n. For V\ G [k, 1] let 5(t>i) = sup{C(f i, z) \ {y i, z) G Aq} [which equals 
C(vi, z(vx)) for some z{yi) such that (v\, z{yi)) G dA^ PI (0, l) 2 since for fixed any v± the function 
it i — y C(v±, u) i increasing] and set 5 = inf^g^i] 5{v i), which is strictly positive due to compactness 



of OAq fl (0, l) 2 and continuity of C. We will now show that this choice of 5 yields (2.3). Now, if 
C(v) < 5, we have either v\ < n [then v G Aq since C(0,t>2) = 0] or v\ > k. In the latter case, 



C(v) < 5(vx) and monotonicity of C imply v G Aq. This proves (2.3) and completes the proof of 
Theorem EH □ 



As a consequence of Theorem 2.2 and the continuous mapping Theorem [see e.g. Theorem 1.3.6 in 



Van der Vaart and Wellner (1996)], we obtain the weak convergence of a corresponding Cramer- 



von-Mises and Kolmogorov-Smirnov type test statistic, i.e. 



T 



n,Lo 



[0,1]' 



{M n (x,y, z)} 2 d(x,y, z) 
T n ,Ks = sup \B. n (x,y,z)\ 

[Q,l]3 



~* T c,u = / {M c (x,y,z)} d(x,y,z) 

J[0,1] 3 

^ Tc jK s = sup \U c {x,y,z)\ , 

[0,1] 3 



(2.4) 
(2.5) 



which can be used to construct an asymptotic test for the hypothesis of associativity. Since 
T n) M - > oo [M G {L 2 , KS}] if the copula is not associative the null hypothesis should be rejected 
for unlikely large values of T nj M- This gives rise to the demand for critical values of Tc,m which 
can be obtained by multiplier bootstrap methods as described in the subsequent paragraph. 
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2.2 A multiplier bootstrap approximation 

It is the purpose of this Section to provide a bootstrap approximation for the distribution of the 
limiting variables Tc,m whose variances depend on the unknown copula in a complicated manner. 
We begin with an approximation of the distribution of the limiting process He*. For this purpose 
we rewrite the decomposition of the process Gc defined in (2.2) as 

Uc{x, y, z) = M c {x, C(y, z)) - C 1 {x, C(y, z))M c (x, 1) - C 2 {x, C(y, z))M c (l, C(y, z)) 

- {n c (C(x, y),z)- C l (C{x, y), z)M c {C(x, y), 1) - C 2 (C(x, y), z)M c {l, z) 

+ C 2 (x, C{y, z)) [M c (y, z) - C 1 (y, z)M c (y, 1) - C 2 (y, z)M c (l, z) 

+ Cx{C(x, y), z) \m c (x, y) - C x {x, y)B c (x, 1) - C 2 (x, y)B c (l, y) } . 

In the following discussion the symbol 



G, 



A G 



(2.6) 
(2.7) 



denotes weak convergence in some metric space D conditionally on the data in probability [see 
Kosorok (2008)]. More precisely, (2.7) holds for random variables G n = G„(X 1 , . . . , X n , £ 1; . . . £ n ), 



G E D if and only if 



and 



sup \E(h(G n )-m(G)\^0 



E^h(G n )* - E^/i(G n )* 4 for every h e BL X 



(2.i 



(2.9) 



where 



{/: 



<l,|/(/3)-/( 7 )|<ci(/3,7)V7,/3eD} 



denotes the set of all Lipschitz-continuous functions bounded by 1. The subscript £ in the ex- 
pectations in (2.8) and (2.9) indicates the conditional expectation with respect to the weights 
£ = (£i,...,£ n ) given the data and h{G n )* and fo(G n )* denote measurable majorants and mi- 
norants with respect to the joint data, including the weights £. Note also that condition (2.8) 
is motivated by the metrization of weak convergence by the bounded Lipschitz-metric, see e.g. 
Theorem 1.12.4 in Van der Vaart and Wellner (1996). 

The process Be can be approximated by multiplier bootstrap methods, see Bucher (2011); Bucher 



and Dette 


(2010 


); 


Remillard and Scaillet 


(2009 


); 


Segers 


(2011 



[ore precisely, let £i, . . . £ n denote 
independent identically distributed random variables with mean and variance 1 such that 



lie 



ill 2,1 



and consider the process 



= / v^G&l > x ) dx < 

oo, 

Jo 



(2.10) 
(2.11) 



7 



where 



n 



denotes a multiplier bootstrap version of the estimator. It was shown in Bucher and Dette (2010) 



and in more detail in Bucher (2011) that 



or -w m>c 



i.e. the process a| defined in (2.11 ) converges weakly to Be in Z°°([0, l] 2 ) conditionally on the data 



in probability in the sense of Kosorok (2008). 



For the approximation of the partial derivatives in (2.6) let C p be some estimator of C p ; for instance 



an estimator based on the differential quotient as in Remillard and Scaillet (2009) defined by 

S(u) := 



62(11) 



C n {u\+h,U2)-C„(u\-h,U2) 
2h 

C n (2h,u 2 ) 
2h 

U2—C„(l—2h,U2) 
2h 

C n Qi ,u 2 +h)-C n (ui ,U2—h) 
2h 

C„(m,2fo) 
2h 

ui— C n («i,l— 2h) 
2h 



if ui G [h, 1 - h] 
if Mi G [0, h) 
if Mi G (1 - h, 1] 

if u 2 G [h, 1 - 
if m 2 G [0, h) 
if w 2 G (1 - h, 1] 



(2.12) 



(2-13) 



where h — h n — >• such that inf„ h n yjn > [for a smooth version of these estimators see Scaillet 



(2005) 



Theorem 2.3. Assume that there exists a constant K such that HC'pHoo < K for alln G N,p = 1, 2 
and that 



sup 

ue[o,i] 2 :M P e[<5,i-<5] 



C p (u) - C p (u) 



-> 



/or a// 5 G (0, 1/2). If moreover the copula C satisfies Condition 2.1 and if the multipliers & satisfy 



(2.10), i/ien £/ie multiplier process El| defined as 



H|(x, y, *) = c4(ar, C ft (y, 2)) - Ci(x, C n {y, z))a{(x u 1) - C-jfo C n (y, z))a{{l, C n (y, z)) 

- {^ n (C n (x,y),z) - Ci(C n (x,y),z)ai(C n (x,y), 1) - C 2 (C n (x, y), *)a£(l, z) 

+ C 2 (x, C n (y, z)) {ai(y, z) - Ct(y, z)ai(y, 1) - C 2 (y, z)a£(l, z) 

+ C 1 (C n (x,y),z) lai(x,y) - C 1 (x,y)al(x, 1) - C 2 (x,y)a&(l,y) 



converges weakly to the process He* conditional on the data in probability, i.e. ~» Mc- 
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Proof. Define the process HT| by substituting the estimators C X ,C 2 and C n in the definition of 
HF| by the true but unknown objects C X ,C 2 and C. By Lemma A.l in Biicher (2011) it suffices to 
show that 



oo = sup \Ml(x,y,z) -B^(x,y,z)\ -> 0. 

(x,j/ 1 z)6[0,l]8 



Using the triangle inequality we have to estimate the following 12 summands 

li - HSIU < ||o4(a:, *)) - a« (x, C(y, z)) |U 

+ ||Ci(x,C n (y,z))c4(x,l)-C'i(x,C , ( 3 / J 2))c4(x,l)|| 00 
+ ||S(£, C n (y, z))a« (1, C n {y, z)) - C 2 (x, C{y, z))a£(l, C(y, z))\\ 

+ \\a&(C n (x,y),z) - C^(C(X,?/),2:)||oc 

+ y), z)ai(C n (x, y),l)- C x (C(x, y), z)a&{C(x, y), 1) || 

+ \\C 2 (C n (x,y),z)ai(l,z)-C 2 (C(x,y),z)ai(l,z)\\ oc 
+ \\C 2 (x,C n (y,z))ai(y,z) - C 2 ( l x,C(y,z))a&(y,z)\\ 00 
+ \\C 2 (x, C n (y, z))C x (y, z)a&(y, 1) - C 2 (x, C(y, z))C x (y, z)af l (y, 1 
+ \\C 2 (x, C n (y, z))C 2 (y, *)a&(l, *) - (%(>, C(y, z))C 2 (y, *)a£(l, * 
+ ||C'i(C n (x,y),z)Q;^(z,y) - C x (C(x,y), z)af l (x,y)\\ 0O 
+ ||Ci(C rt (x, y), z)C x (x, y)afc(x, 1) - C x {C(x, y), z)C x (x, y)a£(x, 1 
+ \\C x {C n {x, y),z)C 2 (x, y)ai(l, y) - C x (C(x, y), z)C 2 {x, y)a|(l, 



of which one of the hardest cases will be considered exemplarily in the following, namely the third 
summand 



sup 

(x,y,z)£[0,l] 3 



C 2 {x, C n (y, z))a£(l, C n (y, z)) - C 2 (x, C(y, z))a£(l, C(y, z)) 



< 



The treatment of the other summands is similar and is omitted for the sake of brevity. We estimate 

C 2 {x, C n (y, z))a£(l, C n (y, z)) - C 2 (x, C(y, «))a£(l, C(y, z)) 

C 2 (x,C n (y, z)) - C 2 (x,C n (y, z)) x |a£(l,C n (y,z))| 

+ C 2 (x,C n (y,z)) -C 2 (x,C(y,z)) x |a£(l, C n (y, z))\ 

+ C 2 (x,C(y,z)) x \ai(l,C n (y,z)) - a&{l,C(y,z))\ 
--:A x (x,y,z) + A 2 (x,y,z) + A 3 (x,y,z) 
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and consider each term separately. For arbitrary e > and S G (0, 1/2) we estimate 



(sup Ax(x,y, z) > e) < P 



sup A\ (x, y, z) > e ) + P 

C n (y,z)E[5,l-S] 



sup Ai(x, y,z) > e 

C n (y,z)f[5,l-S] 



where we suppressed the index (x,y,z) G [0, l] 3 at the suprema. The first probability can be 



made arbitrary sma 
see Theorem 2.3 in 



1 by the assumptions on Ci and by the asymptotic tightness of the process a|, 

Biicher (2011). For the second summand use uniform boundedness of C2 and 
the fact that the (unconditional) limit process B c (l, •) of a|(l, •) is a standard Brownian bridge 
having continuous trajectories which vanish at and 1. By decreasing 5 the probability can be 



made arbitrary small, see Segers (2011) for an rigorous treatment of this argument. 



Since C2 is uniformly continuous if the second coordinate is bounded away from zero and one the 
second summand A 2 (x, y, z) can be treated similarly. Regarding A 3 (x, y, z) note that a| is asymp- 
totically uniformly equicontinuous [Theorem 2.3 in Biicher (2011)] and that sup/^WQ^p \C n (y, z) — 

C(y,z)\ — > which yields 



sup |a£(l, C n (y, z)) - a£(l, C(y, z)) \ -> 0. 

fe^)6[0,l] 2 



F 

By boundedness of C2 this yields the assertion sup^^wo^p A 3 (x,y,z) — > 0. 



□ 



Remark 2.4. 

(a) Note that the assumptions on the estimator C p for the partial derivatives C p are e.g. satisfied 
for the estimators defined in (2.12) and (2.13), see Lemma 4.1 in Segers (2010). 



(b) Note that Theorem 2.3 holds independently of the hypothesis of associativity. As a consequence 
of the continuous mapping theorem for the bootstrap, see Proposition 10.7 in Kosorok (2008), we 
can conclude that 



TP? 



{M.i(x,y,z)} d(x,y,z) ~» T CjL2 , 
[o,i] 3 « 



T 



n,KS 



sup \mi{x,y,z)\ -w T CiKS (2.14) 



[0,1] 3 



and the latter convergence suggests to use the following approach in order to obtain an asymptotic 
level-a test for the hypothesis of associativity. 

1. Compute the statistic T niM [M G {L 2 ,KS}]. 

2. Choose the number of bootstrap replications B G N. For b = 1, . . . , B simulate independent 
replications of the random variables £1, . . . ,£„ and denote the result form the 6-th iteration 

by £l,6, • • • ,£ n ,&- 



3. For b = 1, .. . ,B compute the statistics T^u defined in (|2.14[) from the data X 



. Al- 



and the multipliers £1^, . . . , £ n> & and determine the (1 — a)-quantile q{_ aM of the empirical 



distribution of the sample {T^J^} 



b=l,...,B- 
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4. Reject the null hypothesis of associativity whenever T n _M > q{_ a M 
p t p 

Since T n) M — > oo and T„ M ~> Tc,m if the copula is not associative the test is consistent against 



all alternatives satisfying the conditions of Theorem 2.3 



3 Testing Archimedeanity 

As stated in the Introduction a bivariate copula C is Archimedean if and only if C is an associative 
copula satisfying C(u,u) < u for all u G (0, 1). Associativity has been dealt with in the preceding 
paragraph and it remains to handle non- Archimedean copulas which my be associative but satisfy 



C(q, q) = q for some q G (0, 1). Due to Theorem 1 in Jaworski (2010) or by the results in Section 



2.4 of Alsina et al. (2006) all those copulas may be expressed as an ordinal sum of Archimedean 



copulas. An ordinal sum copula is defined as following [cf. Section 3.2.2 in Nelsen (2006)]: let 
{Ji}iei be a countable partition of non-overlapping closed intervals Jj = [a*, &,] whose union is 
[0, 1]. If moreover {Cj} ie / is a collection of copulas, then the ordinal sum of {Ci} ie i with respect 
to { Ji}i<=i is the copula C defined by 



C(u) 



min{M 1 ,M 2 } otherwise. 



Note that C(bi, hi) = b{ for all 6, and that ordinal sum copulas put no mass on [0, l] 2 \ |J igJ Jj x Jj. 
In Figure [T] we illustrate the ordinal sum of a Gumbel copula C\ with parameter Q\ = 1.5 and a 
Clayton Copula C% with parameter 62 = 1, where J\ = [0, 1/2], J2 = [1/2, 1]. Note that Kendall's 
r of C is equal to 2/3, while it equals 1/3 for both C\ and C 2 . 

In order to check for C(q, q) = q for some g G (0,1) we propose the following modification of the 
statistic T Ut M 

§n,M = TT n>M + k n (j)(A n (C n )), 



where T nj M is defined in (2.4) and (2.5), k n ~ n a ,a G (0,1/2) is some constant chosen by the 



statistician, is some increasing function with 0(0) = and 

A n (C n ) = m^\-(l--) : cJ -,-)=- 

Intuitively, A n {C n ) should be "large" for copulas which satisfy C(q, q) = q for some q G (0, 1). For 
a decent choice of k n and we refer the reader to Section |4| 

In Figure [2] we illustrate the points i/n at which C n (i/n,i/n) = i/n for two specific examples, 
the Clayton copula with 9 = 1 and the ordinal sum copula depicted in Figure [Tj The solid and 
dashed lines correspond to the true copula and the empirical copula calculated for a set of n = 100 
observations, respectively. For the ordinal sum copula there always exist some points i/n in a 
neighbourhood of 1/2 such that C n (i/n,i/n) = i/n, see the proof of the following Proposition, 
which is sufficient for the derivation of the asymptotic properties of the statistic S^m- 
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Copula 



Simulated Observations 




Figure 1: Left picture: Ordinal sum copula. Right picture: 500 corresponding simulated observa- 
tions. 

Proposition 3.1. 



a) Suppose C is an Archimedean copula satisfying Condition 2.1 and that the coefficients of tail 
dependence 

C(u,u) l-2u + C(u,u) 

al = lim and Ajj = Iim 



w->o u u-+i 1 — u 

exist and are smaller than 1. Then it holds 

A n {C n ) = o P (n- Q ) 

for any a < 1/2. 
b) If there exists a q G (0, 1) such that C(q, q) = q then it holds 

A n {C n ) > g(l-g) + op(l). 

Proof. We start with the proof of a). First choose 5 > and A < 1 such that 

C(u,u) l-2u + C(u,u) 

V ^ A 

u 1 — u 

for all u G [0, 5] U [1 — 8, 1] and use the decomposition 

A n (C n ) = A n (C n , [0, 5)) + A n (C n , [5, 1 - 5}) + A n (C n , (1 - 5, 1]), (3.1) 

12 



Clayton-Copula (tau=1/3) 



Ordinal-Sum-Copula (tau=2/3) 




Figure 2: The solid lines show the diagonal section of a Clayton copula (left) and an ordinal sum 
copula (right), while the dashed line show one realization of the corresponding empirical copula. 
The circled points mark the locations where C n {i/n,i/n) = i/n. 

where 

A n (C n , B) = max {i/n (1 — i/n) : C n (i/n, i/n) — i/n and i/n G B} 

for some set B C [0, 1] (with the convention that max0 = 0). Consider each term separately and 
define M n = sup u6 r 0)1 ] \C n (u,u) — C(u,u)\ , which is of order Of(n~ l l 2 ) under Condition 
let i/n G (0,5) be such that C n (i/n,i/n) = i/n. Due to the estimate 



2.1 



Now 



z'/n(l - A) < i/n ( 1 - C<yl l n ' ) = i/ n - C(i/n,i/n) = C n (i/n,i/n) - C(i/n,i/n) < M n 



i/n 

we have i/n(l — i/n) < i/n < M n /(1 — A) and we can conclude that 

A n (C n , [0,5)) < ^ = P (n" 1/2 ). (3.2) 

A similar calculation shows that for i/n G (1 — 5, 1] with C n {i/n, i/n) = i/n we have (1 — i/n)(l — 
A) < M n which in turn implies 

A n {C n , (1 - 5, 1]) < ^ = P (n-^ 2 ). (3.3) 



It remains to estimate the second summand A n (C n , [5,1 — 5]) of the decomposition (3.1). For 



continuity reasons we can choose a k > such that u — C(u,u) > k for all u G [5, 1 — 5]. If there 
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was a q G [5, 1 — 5} such that C n (q, q) = q, it would follow that M n > C n (q, q) — C(q,q) > k and 
therefore we have for any e > 

P (n a A n (C n , [5, l-S\)>e)< P(3 g G [5, 1 - 6} : C n (g, g) = g) < P(M n > k) 4 0. 



A combination of (3.2) and (3.3) with this result proves part a) of the proposition. 
For the proof of part b) let m = #{1 < i < n : (F^X^), F 2 (X i2 )) G [0,g] 2 }. Since C(q,q) = q 
implies that the mass of C is concentrated on (0,g) 2 U (g, l) 2 we have (Fi(Xn) , F 2 (X i2 )) G [0,g] 2 
if and only if < X ni - n ^ and X i2 < X ni:n ^, where Xj :ntP = F~ p (j/n) denotes the j-th order 
statistic of X\ p , ...,X np (for p = 1,2). This yields C n (ni/n,ni/n) = n\/n, which entails the 
assertion by 



A n {C n ) > — [ I 
n 



71 



-> q(l - q) > 0. 



□ 



Remark 3.2. 



a) The conditions on the coefficients of tail dependence in part a) of Proposition 3.1 can be 
equivalently expressed by conditions on the regular variation of the Archimedean generator of 



C. For a thorough discussion of these issues the reader is referred to the work of |Charpentier 
|and Seg" 



?crs 



(2009). 



b) Exploiting Theorem G.l in Genest and Segers (2009) and Proposition 4.2 in Segers (2011) 



one can improve the rate of convergence in part a) of Proposition 3.1 to any a < 3/4. It is 
our conjecture that the term is in fact of order Op(l/n), but we were not able to derive the 
asymptotic distribution of nA n {C n ). Since we do not need a refined rate for our purposes here 
we omit a deeper discussion and defer these issues to future research. 



From now on, suppose that the conditions of Theorem 2.3 and Proposition 3.1 hold. We can 
conclude that § ni M weakly converges to Tc,u if t ne copula C is Archimedean, while § nj M converges 
to +oo in probability if C is non-Archimedean, i.e. if it is either non-associative (by the results of 



Section [2]) or if there exists a q G (0,1) such that C(q,q) = q (by Proposition 3.1). The quantiles 
of Tcm can be approximated by the multiplier method described in Section 2.2 Analogously 



to the discussion at the end of the Section 272] we can use the multiplier bootstrap to obtain 
an asymptotic level-a test for the hypothesis of Archimedeanity, which is consistent against all 
alternatives satisfying the Condition |2.1 Its finite sample properties will be investigated in the 
following section. 



4 Finite sample properties 

We conclude this paper with a simulation study regarding the finite sample performance of the 
proposed tests for Archimedeanity and Associativity. We consider the following six copula models: 
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The Gumbel copula, which is Archimedean. 
The Clayton copula, which is Archimedean. 

The t-copula with fixed degree of freedom df = 1, which is non- associative. 



• The asymmetric negative logistic model [see Joe (1990)] with fixed parameters ip\ = 2/3, ip2 = 
1, which is non-associative. 

• An ordinal sum model based on the partition J\ = [0, 1/2], J 2 = [1/2, 1] and the Gumbel 
(Ci) and Clayton (C 2 ) copula, denoted by Ordinal^. The model is associative. 

• An ordinal sum model based on the partition J\ = [0, 1/2], J 2 = [1/2, 1] and the two Clayton 
(Ci = C 2 ) copulas, denoted by Ordinal^. The model is associative. 

The parameters of the models are chosen in such a way that the coefficient of upper tail dependence 
\u is either 1/4 or 1/2 [for the asymmetric negative logistic model] or that Kendall's-r is either 
1/3 or 2/3 [for the remaining five models]. For ro r dinai A = 1/3 (or 2/3, resp.) we chose rcumbci = 
(1/3) and r C i ayton = -2/3 (1/3), while r clayton = -1/3 (1/3) for r 0r dinai s = 1/3 (2/3). 
We generated 1000 random samples of sample sizes n = 200 and n = 500 and calculated the 
empirical probability of rejecting the null hypotheses of Archimedeanity or Associativity for M e 
{L 2 ,KS}. For each sample of size n = 200,500 we carried out B = 200 Bootstrap replications 
based on the multiplier method, where we chose a U({0, 2})-distribution for the multipliers [i.e. 
P(f = 0) = P(f = 2) = 1/2, s.t. fi = r = 1] and used h = n~ 1/4 to estimate the partial 
derivatives. The critical values of the tests are determined by the method described in Section 
2] The penalty term S n) M — T n) M = k n (p(A n (C n )) is chosen in the following data-adaptive way: 
first of all, we set <p(x) = (Ax) 2 in order to give more emphasis to values around the maximal 
value of A n (C n ) [which equals 1/4]. The constant k n is chosen according to the distribution of the 
bootstrap approximation: if <7oo5M denotes the 0.05-quantile of the sample {^ u }b=i,...,B we set 
05 M n • T ne l^ter choice guarantees that under H the error term is small compared to 
the distribution of Tq,m- 

The results are stated in Table [TJ The entries of the table represent the empirical probabilities 
of rejecting the null hypothesis of Archimedeanity and of Associativity [in brackets] for both the 
L 2 -test [first two columns] and the KS-test [last two columns]. We observe that the nominal level of 
the four tests are accurately approximated for the four Archimedean copulas under investigation. 
The L 2 -test tends to be more conservative than the KS-test. Also note that the values for § ni M 
and T n _M differ only by a very small amount meaning that the penalty term k n (p(A n (C n )) is of 
negligible magnitude under the null hypothesis. 

The t-copula and the asymmetric negative logistic models are non-associative and the results in 
Table [T] reveal that these deviations are detected by both tests for Associativity, with better results 
for the i-copula and for stronger dependence [measured by either r or Xjj]. The power properties 
of the L 2 -test outclass the properties of the KS-test for all four non-associative models under 
investigation, such that the former test seems to be generally preferable. 
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Regarding the (associative) ordinal sum models both tests for associativity are very conservative. 
Note that the asymptotic theory developed in Section [2] does not apply for these models since the 
partial derivatives of the corresponding copulas are not continuous on {1/2} x [0, 1] and [0, 1] x 
{1/2}. Regarding the power properties the KS-test for Archimedeanity performs slightly better 
for the ordinal sum alternatives. 

Acknowledgements This work has been supported by the Collaborative Research Center "Sta- 
tistical modeling of nonlinear dynamic processes" (SFB 823, Teilprojekt Al, CI) of the German 
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L 2-Test KS-Test 
0.1 0.05 0.1 0.05 



n=200: 


Clayton(r = 


1/3) 


0.071 


(0.071) 


0.038 


(0.037) 


0.088 


(0.088) 


0.050 


(0.050) 


Glayton(r = 


2/3) 


0.030 


(0.016) 


0.011 


(0.009) 


0.124 


(0.078) 


0.068 


(0.036) 


G umbel (r = 


1/3) 


0.079 


(0.079) 


0.043 


(0.043) 


0.082 


(0.082) 


0.046 


(0.045) 


Gumbel(r = 


2/3) 


0.034 


(0.032) 


0.015 


(0.013) 


0.108 


(0.098) 


0.065 


(0.057) 


t(r= 1/3, df 


= 1) 


0.953 


(0.953) 


0.886 


(0.884) 


0.562 


(0.558) 


0.380 


(0.376) 


t(r = 2/3, df 


= 1) 


0.748 


(0.726) 


0.592 


(0.564) 


0.392 


(0.355) 


0.258 


(0.231) 


Aneglog(A l/ = 


= 0.25) 


0.112 


(0.112) 


0.061 


(0.061) 


0.105 


(0.105) 


0.059 


(0.059) 


Aneglog(At/ = 


= 0.5) 


0.641 


(0.641) 


0.536 


(0.536) 


0.363 


(0.356) 


0.225 


(0.222) 


Ordinal^ (r = 


= 1/3) 


0.996 


(0.000) 


0.827 


(0.000) 


1.000 


(0.012) 


1.000 


(0.005) 


Ordinal a (t = 


= 2/3) 


1.000 


(0.004) 


1.000 


(0.001) 


1.000 


(0.079) 


1.000 


(0.041) 


Ordinal^ (t = 


= 1/3) 


1.000 


(0.000) 


1.000 


(0.000) 


1.000 


(0.045) 


1.000 


(0.021) 


Ordinal B (r = 


= 2/3) 


1.000 


(0.006) 


1.000 


(0.001) 


1.000 


(0.057) 


1.000 


(0.030) 


n=500: 


Clayton(r = 


1/3) 


0.082 


(0.082) 


0.051 


(0.051) 


0.088 


(0.088) 


0.036 


(0.036) 


Clayton(r = 


2/3) 


0.062 


(0.059) 


0.032 


(0.027) 


0.090 


(0.082) 


0.046 


(0.039) 


Gumbel(r = 


1/3) 


0.091 


(0.091) 


0.045 


(0.045) 


0.105 


(0.105) 


0.050 


(0.050) 


Gumbel(r = 


2/3) 


0.072 


(0.070) 


0.033 


(0.032) 


0.124 


(0.121) 


0.068 


(0.066) 


t(r =1/3, df 


= 1) 


1.000 


(1.000) 


1.000 


(1.000) 


0.954 


(0.953) 


0.871 


(0.871) 


t(r = 2/3, df 


= 1) 


0.998 


(0.990) 


0.998 


(0.990) 


0.818 


(0.811) 


0.655 


(0.650) 


Aneglog(A l/ -- 


= 0.25) 


0.237 


(0.237) 


0.173 


(0.173) 


0.124 


(0.124) 


0.069 


(0.069) 


Aneglog(At/ = 


= 0.5) 


0.979 


(0.979) 


0.947 


(0.947) 


0.716 


(0.716) 


0.584 


(0.584) 


OrdinaU(T = 


= 1/3) 


1.000 


(0.000) 


0.980 


(0.000) 


1.000 


(0.022) 


1.000 


(0.007) 


OrdinaU(T = 


= 2/3) 


1.000 


(0.021) 


1.000 


(0.009) 


1.000 


(0.093) 


1.000 


(0.038) 


Ordinal^ (r = 


= 1/3) 


1.000 


(0.000) 


1.000 


(0.000) 


1.000 


(0.005) 


1.000 


(0.023) 


Ordinals (r = 


= 2/3) 


1.000 


(0.013) 


1.000 


(0.004) 


1.000 


(0.082) 


1.000 


(0.037) 



Table 1: Simulated rejection probabilities for the multiplier bootstrap-tests for Archimedeanity and 
for Associativity (in brackets). The sample size is n = 200 or n = 500, B = 200 Bootstrap- 
replicates and 1000 simulation runs have been performed. The first four lines are Archimedean 
copula models, the t-models are not associative and the ordinal sum-models are associative, but not 
Archimedean. 
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